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SUMMARY 


Nonlinear equations of motion of symmetrical ly laminated anisotropic 
plates are derived accounting for von Karman strains. The effect of trans- 
verse shear is included in the formulation and the rotatory inertia effect 
is neglected. Using a single-mode Galerkin procedure the nonlinear modal 
equation is obtained. The direct equivalent linearization method is em- 
ployed for solution of this equation. The response to acoustic excitation 
of moderately thick composite panels is studied. Further, the effects of 
transverse shear on large deflection vibration of laminates under random 
excitation are studied. Mean-square deflections and mean-square inplane 
stresses are obtained for some synmetric graphite-epoxy laninates. Using 
equilibrium equations, and the continuity requirements, the mean-square 
transverse shear stresses are calculated. The results obtained will be 
useful in the sonic fatigue design of composite aircraft panels. The analy- 
sis is presented in detail for simply supported plates. The analogous equa- 
tions for a clamped case are given in the appendix. 

INTRODUCTION 

Acoustically induced fatigue failures in aircraft structures have been 
a design consideration for the past three decades. With the advent of the 
jet engine which produced high intensity acoustic pressure fluctuations on 
aircraft surfaces, the problem acquired prominence. The number of acoustic 
fatigue failures have resulted in unacceptable maintenance and inspection 
burdens associated with the operation of aircraft. Therefore, accurate 
design methods are needed to determine the acoustic fatigue life of struc- 
tures. Numerous analytical studies Refs. 1-13 and experimental investiga- 
tions Refs. 7, 14-22 on sonic fatigue design of aircraft structures have 



been undertaken during the past decade to help in providing the needed in- 
formation. 

The majority of analytical studies on flat panels to date have been 
formulated within the framework of linear or classical plate theory which 
assumes small deflections. Current analytical design methods Refs. 7, 9, 

12, 13 for sonic fatigue prevention are based essentially on linear struc- 
tural theory. Test results on various aircraft panels reported in the 
literature. Refs. 7, 13-17, 19-22 however, have shown that high noise levels 
produce nonlinear large deflection behavior in such panels. Recently, 
analytical efforts Refs. 1-3, 5, 6, 8, 16 have demonstrated that the 
prediction of panel random response is greatly improved by including the 
large deflection effects in the formulation. In all these efforts, both 
analytical and experimental, the thickness of panel is very small as 
compared with its length (a/h > 100). The effect of acoustic excitation on 
moderately thick composite structural panels has not been investigated. 
Further, the effects of transverse shear on large deflection vibrations of 
laminates under random excitation have also not been studied. 


The classical theory of plates, in which it is assumed that normals to 
the midplane before deformation remain straight and normal to the plane 
after deformation, under predicts deflections and over predicts natural 
frequencies and buckling loads. Such results are due to the neglect of 
transverse shear strains in the classical plate theory. The errors in de- 
flections, stresses, natural frequencies and ouckling loads are even higher 
for plates made of advanced composites like graphite-epoxy and boron-epoxy, 
whose elastic modulus to shear modulus ratios are very large (e.g., of the 
order of 25 to 40, instead of 2.6 for typical isotropic materials). These 
high ratios render classical theories inadequate for the analysis of compos- 
ite plates. Many plate theories exist that account for transverse shear 
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strains Refs. 23-34. Recently, the dynamic von Karman plate theory has been 
extended to include the shear and rotatory inertia effects by Sathyamoorthy 
and Chia for nonlinear free vibrations of anisotropic rectangular (Ref. 35) 
and Skew (Ref. 36) plates. Sivakumaran and Chia have extended this approach 
to generally laminated anisotropic thick plates. Ref. 37. 

In this report, the equations of motion are derived from plate theory 
which takes Von Karman large deflection strain-displacement relations into 
account. The transverse shear deformation effects are included and the 
rotatory inertia effects are neglected. The system of equations is then 
simplified to two coupled nonlinear differential equations in terms of 
transverse displacement and a stress function. Due to the complex nature of 
the problem, the study is restricted to a single-mode response. A deflec- 
tion function that represents the first mode is assumed; and corresponding 
to the assumed mode, a stress function satisfying the different inplane edge 
boundary conditions is obtained by solving the compatibility equation. The 
Galerkin method is applied to the governing equation of motion in the trans- 
verse direction using the assumed displacement function as the weighting 
function. This yields a nonlinear, nonhomogeneous, second-order differen- 
tial equation of the response in time. Slope functions that include the 
transverse shear deformation effects are assumed and the coefficients are 
evaluated by Galerkins approximation. Finally, the random responses from 
cases based on four formulations at various acoustic loadings are evaluated 
for simply supported rectangular symmetrical laminates. The four cases of 
formulations are: 1) linear, small deflection plate theory without trans- 
verse shear deformation, 2) linear theory with shear, 3) large deflection 
without shear and 4) large deflection with shear. The excitation is assumed 
to be stationary, ergodic and Gaussian with zero mean; its magnitude and 
phase are uniform over the plate surface. 
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The equivalent linearization method is employed. Root-mean-square (RMS) 
deflections, RMS inplane stresses and frequencies are calculated. Using 
three-dimensional equilibrium equations and the continuity requirements, the 
RMS transverse shear stresses in the laminate are estimated. 

EQUATIONS OF MOTION 

Consider an initially flat, rectangular, elastic plate of constant 
tnickness h in the z-direction, length a in the x-direction and width b 
in the y-direction, see Fig. 1. The reference plane z=0 is located at the 
undeformed middle plane. 

The displacement components that include the effects of transverse 
shear deformation are assumed in the form (Ref. 38): 

U(x,y,z,t) = u(x,y,t) + z ct(x,y,t) 

V(x,y,z,t) = v(x,y,t) + z B(x,y,t) (1) 

W(x,y,z,t) = w(x,y,t) 

in which U, V and W are the inplane and transverse displacements in x,y 
and z directions, respectively, u, v and w are the values of U, V, W 
at the midsurface of the reference plane, and a, 6 are slope functions in 
the xz and y z planes due to bending only. These are averaged components 
of direction change of the normal to the undeformed middle surface. The 
total strains for the laminated plate can be expressed as 


4 


e 

x 


e 


y 


e 


xy 


e + z a, 
x x 


e y + ze *y 


e xy + z(a ’y + 6 >x ) 


( 2 ) 


e 

z 


= 0 


e 


xz 


0 + w, 

X 


e 


yz 


e + w. 


y 


where e , e and e are the engineering shear strains, 
xy xz yz 

The von Karman large deflection strain-displacement relations are given 
by 
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As in the classical plate theory, stress resultants and stress couples are 
defined as 
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where Q's are the transformed reduced stiffness. Thus for symmetrical 
1 aminates 
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and elements of laminate stiffnesses A., and D.. 

ij 


are defined as 
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Partial inversion of equation (6) and inverted form of equation (7) can be 
written as 
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The plate equations are obtained by considering the equilibrium of an 
element of the kth layer of the laminate. Integrating the equilibrium 
relations over the plate thickness h, neglecting inplane inertia terms and 
retaining the nonlinear terms in accordance with the Von Karman assumptions, 
leads to the following equations of motion 


N + N = 0 
x,x xy,y 


N + N = 0 
xy,x y,y 


M + M = q 
x,x xy,y ^x 


M + M = Q 
xy,x y,y v y 


< N x w ’x + N xy w, y^’x + ^ N xy w, x + N y W *y^’y 


+ l + 1 „ + P(t) = p h w 

A i A JlJ 


Expanding equation (11) and introducing a tracing constant T , we can 


a + w, x = T s (A 55 * Q x + A 45 * Q v ) 
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8 + W ’y " T $ (\ 5 * \ + \ 4 * Qy) 


( 18 ) 


in which T $ takes the value of either 1 or 0. If transverse shear 
deformation effects are neglected T $ = 0 and equations (17) and (18) 
reduce to 

a = -w, x and B - - w, y (19) 

Using equations (14) and (15) and making use of constitutive relations given 
in equations (10), equations (17) and (18) can be shown to be 

a + w, x = bia, xx + b 2 a, xy + b 3 a, yy + M, xx + b 5 (3, xy + b 6 8, yy (20) 

6 + w, y = b 7 a, xx + b 8 a, xy + bga, yy + bi 0 8, xx + bnB, xy + bi 2 8, yy (21) 

where the coefficients b i are defined as 

bi = T $ (A55* On + A45* Dig) 

b2 = T s (2 A55* Di 6 + A45* (D12 + D 66 ) ) 

b 3 = T s (A55* D 6 6 + A45* D26) 

bi, = T $ {A55* Di 6 + A45* D6 6 ) 

bs = T s (A55* (D12 + De 6 ) + 2 A45* D2 6 ) 
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b 6 - T s (Ass* »26 ♦ * 45 * D 2 2 ) 


b7 = T s (A 45 * D 11 + A 44 * Di 6 ) 

bs = T s ( 2 A 45 * Dl6 + A 44 * (D 12 + D66)) (22) 


bg = T- (A 45 * D66 + A 44 * D 26 ) 


bio = T 5 (A45* Dig + A44* Dg g ) 


b ll = T s (*45* (°12 + D 6e) + 2 A 44 * D 2 6 ) 


bl2 = T s (A45* D 26 + A 44* D 22 ) 


Equations (20) and (21) can be written as 


w, x + J(a) + K(B) = 0 


(23) 


w, y + L(a) + M(e) = 0 


(24) 


where J,K,L and M are the operators defined as 
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Solving equations (23) and (24) we get 


M(w, x ) - K(w, y ) = N(a) 


J(w, y ) - L(w, x ) = N(6) 


(25) 

(26) 


where the operator N = KL - MJ can be shown to be 
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with coefficients k.. defined as 


ki = b4 b7 - bi bio 


k 2 = b»» bs + bs b7 - bi bn - b2 bio 


k3 * bi* bg + bs bs + bg b7 - bi bx 2 - b2 bn - b3 bio 


11 



( 28 ) 


k 4 b 5 b 9 + b 6 b 8 - b3 bll - b2 bl 2 
k5 = be bg - b3 bl 2 

k6 = bi + bio 

k 7 = b 2 + bn 

kg = b 3 + bn 

Using equations (14) and (15) and the relationships given in equations (6) 
and (8) we can get 
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The Airy's stress function F is defined such that 


[N] T = f F »yy F, XX -f. Xy J 


(31) 
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using equation (31) and the relationships given in equations (6) and (8), we 
can get 


N x W -XX + 


N w, + 

y yy 


2N xy w *xy ‘ ♦< F ’ W > 


(32) 


Substituting equations (29) and (32), equation of transverse motion as given 
in equation (16) can be rearranged to be 

Ii + I 2 + I 3 * 0 (33) 


where 


11 = p(t) - p h w 

1 2 = ♦(F.w) (34) 

1 3 = Lj io,B) 

By making use of the operators defined in equations (25) and (26), a and B 
can be eliminated from equation (33) and can be written as 

N(Ii + I 2 ) + Uj (w) = 0 (35) 

where N is the operator as defined in equation (27) and 
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where P^ are defined as 


Pi = b 7 Di6 - b 10 Du 


P 2 = (b 4 - bn) Du + ( b 8 - bi - 3bio) Di6 + b 7 (Di2 + 2D6 6) 


P 3 = ( b 5 - b 12 ) Du + ( b 9 ■ b 2 + 3b 4 - 3bi 1 ) Di6 


+ (bs - bi - bio) (D 12 + 2D66) + 3b7 D 26 


p 4 - b 6 D ii + ( 3bs _ 3bi 2 - b3) Di6 


+ (b 4 - bn + bg - b2) ( Di 2 + 2D6 6) 
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P 5 = 3 b6 Di6 + (bs - bi 2 - b 3 ) ( Di 2 + 2 D6 6) 


14 


+ (bit " ^11 + ^bg ■ 3b 2 ) D 2 6 + (b 8 - bj) D 2 2 


Pe = b 6 (D12 + 2066 ) + ( bs - bi2 - 3 b 3 ) D26 + ( bg - b2) D22 


P 7 = bg D 26 - bg D 22 

Thus the equation of motion in the transverse direction as shown in 
Eq. (35) can be expressed in terms of out-of-plane deflection w, the 

stress function F and the operators N and Ui. Note that with tracing 
constant Ts = 0, the b^ , k^ , and P i are zero, then the operators N 
and Uj reduce to 


N = - 1 
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The compatiblity equation is derived from equation (3) and can be written 
as 
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Maxing use of equations (10) and (31), one can write Eq. (38) as 

A 22 * F, xxxx ' 2A26 * F >xxxy + ( 2A 12 * + A 66*)F, XX yy - 2Ai6*F, xy yy + 

A 1 1 F »yyyy xy ^*XX ^*yy (39) 

Equations (35) and (39) are the governing equations which will be solved by 
employing Galerkin's approach and equivalent linearization method. 

DEVELOPMENT OF SOLUTIONS 

Modal Equation 

Consider a simply supportedt, rectangular, symmetric composite plate of 
dimensions a x b x h with the origin located at the center of the plate. 
The out-of-plane boundary conditions'! are 

x = ± a/2 : w = M x = B = 0 (40) 

y = ± b/2 : w = M^ = a = 0 

For the inplane condition of zero shear stresses at the edges, the 
deflection function is assumed as 

w(x,y,t) = q(t)h cos — cos — (41) 

a b 

The slope functions a and B are assumed as 

tTne analogous equations for clamped case are shown in appendix. 
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( 42 ) 


a(x,y,t) = B! q(t) sin — cos — 

a b 


8(x,y,t) = B 2 q(t) cos ™ sin 

a b 


(43) 


where the constants Bi and B 2 will be evaluated later. Substituting 
equation (41) in equation (39) and solving it, the stress funciton F is 
obtained as 


F = F + F 
c p 


(44) 


in which the particular solution is given as (Ref. 6) 


f. -£!£! (f, . ,00**2) 

H 32 a b 


(45) 


where 
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'10 F 0 1 = ; r - a/ b 

A *22 A n *r 4 


(46) 


The complementary solution F c will now be obtained such that it satisfies 
inplane boundary conditions. For movable edges, the inplane boundary 
conditions are 


b/2 


x = ±a /2 : F >xy-/-b/2 F ’yy = 0 


(47) 
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y = ±b/2 


F, yy = ' 


a/2 

-a/2 


F ’xx dx = 0 


By using the above conditions, it can be shown that F c is zero for movable 
inplane edges. For immovable edges, the inplane boundary conditions of zero 
shear stresses and zero normal displacement at the four edges are 


x = ±a/2 : F, xy = f / (e° - 1 w, x ) dxdy = 0 


(48) 


y = ±b/2 : F, xy = f / (e y - w , y ) dxdy 


The complementary solution is assumed as 


F 


c 



xy 


(49) 


Upon using Eqs. (10) and (41) and enforcing the conditions of Eq. (48), N , 

N and N in Eq. (49) are obtained as 
y xy 


I = ( * 22 k - hi) 
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i = q 2h2ir2 ( h±*~ hi* \ 

y 8(A n *A 22 * . A 12 * 2 ) \ b 2 a 2 / 


(50) 
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The particular solution F p has been obtained and given in equation (45). 
The total stress function, therefore, is F = F p + F c for immovable inplane 
edge conditions. With the assumed deflection w given by equation (41) and 
the stress function F given by Eq. (44), Eq. (35) is then satisfied by 
applying a modified Gal er kin's method: 


/ / [N (Ii + I 2 ) + Ui ( w) ] wdxdy 


+ f% x = -a/ 2 (a) x = -a/2 dy 
b/2 

+ -b/2 ^ x=a/2 ^ x=a/2 ^ 
a/2 

+ ; -a/2 ( V y =- b/2 (6) y - b/2 ^ 


a/2 

+ J -a/2 ( V (B) y^ 2 dX = ° 


which yields a modal equation of the form 


q + a)2 q + x q 3 = 
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(51) 


(52) 
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and m is the modal mass, it is given by 


0 ' • k *£) 4 + ‘•ft 
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Note that with T = 0, k. and P. are zero and Eqs. (53) and (57) 

S I 1 

reduce to 
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ph 2 ir 2 
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where u> is linear radian frequency, X is nonlinearity coefficient, X 
op c 

is an addition to the nonl inearity coefficient due to immovable inplane edge 
conditions and m is the modal mass. For movable inplane edge conditions 
X c = 0. For linear small deflection theory, X = 0. Equation (52) repre- 
sents a single-mode, undamped, large amplitude modal equation with trans- 
verse shear effects taken into account (T s =l). This nonlinear modal equa- 
tion will be solved by employing the method of equivalent linearization. 
Random Response 

It is known that damping has significant effects on the response of 
structures. Therefore, the precise determination of the damping coefficient 
of a structure should be emphasized. The values of damping ratio 
c (=c /c cr ) generally range from 0.005 to 0.05 for the common type of 
composite panel used in aircraft construction (Refs. 12, 20, 22). Once the 
damping ratio is determined from experiment or from existing data, Eq. (52) 
can be expressed in a general form as 

2 3 p ( t ) 

q + 2? u> o q + w o q + Xq = Z±-L 


m 


(60) 



The method of equivalent linearization is used to obtain an approximate 
mean-square amplitude of Eq. (60). 

The basic idea of the equivalent linearization method (Refs. 39, 40) is 
that the approximate response can be obtained from the linearized equation 

q + 2 c<o q + ft q = (61) 

m 

where ft is an equivalent linear or nonlinear frequency. The error of 
linearization, a random process, is 

« = (“o - n2 )d + Ac ! 3 (62) 

which is simply the difference between Eq. (60) and Eq. (61). The method of 

attack is to minimize this error in a suitable way. The classical choice is 

2 

to minimize the mean-square error Ef6 ], that is 

— ~ Ef^ 2 ] = 0 (63) 

9(S! ) 

If the acoustic pressure excitation p(t) is stationary Gaussian, is 
ergodic, and has a zero mean, then the approximate displacement q computed 
from the linearized Eq. (61), is also Gaussian and approaches stationarity; 
this result is due to the fact that the panel motion is stable. Substitut- 
ing Eq. (62) into Eq. (63) and interchanging the order of differentiation 
and expectation yields 
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(64) 


(«* - Q 2 ) E [ q 2] + X E[ q 4 ] = 0 


which leads to the relatonship between the equivalent linear frequency and 
the mean-square displacement as 


2 

n 


2 r 2 ! 

u + 3X Ei q 1 


(65) 


where E[q ] is the mean-square maximum deflection of the laminated compos- 
ite plate. 

The mean-square response of the modal amplitude from Eq. (61) is 


E[q ] = /q S(u>) |H(u>) | do> 


( 66 ) 


where S(w) is the pressure spectral density (PSD) function of the 
excitation p(t). The frequency response function H(u>) is given by 



m(fi - a + ZiCtu jd) 


(67) 


For lightly damped (? < .05) structures, the frequency response curve will 
be highly peaked at the equivalent linear frequency £2 (not at w o as in 
the small deflection linear theory). Integration of Eq. (66) can be greatly 
simplified when the spectral density of the excitation is slowly varying in 
the neighborhood of P, and S(fl.) can be treated as constant in the fre- 
quency band surrounding this nonlinear resonance peak fi; then Eq. (66) 
yields 
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TrS(fi) 


( 68 ) 


Eh 2 ] = 

2 2 
4m ?oj 

o 

To convert the PSD function from (Pa) 2 /radian to (Pa) 2 / Hz, substitute 

n = 2nf 

S(B) = HH (69) 

2ir 

into Eq. (68); then the mean-square deflection becomes 

ET q 2 l = (70) 

8m 2 ?o) n 2 

o 

2 2 

The PSD function S(f) has the units (Pa) /Hz or (psi) /Hz. 

Slope Functions 

In deriving the governing equations of motion, the slope functions a 
and 8 were eliminated as such, and for the determination of linear and 
nonlinear frequencies and the mean-square displ acenent, the slope functions 
need not be known. But for the determination of strains and stresses, slope 
functions a and 8 are to be known. The boundary conditions for a simply 
supported plate are given by Eq. (40). The slope functions given in Eqs. 
(42) and (43) satisfy the boundary conditions. The constants Bj and B 2 
are determined by applying Galerkin method. Using Eqs. (25), (26), and (41) 
and with slope functions a and 8 as the weighting functions, one can 
get: 
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( 71 ) 


jD/2 j-a/2 [ N { a ) _ M ^ j + K ( w )] a dxdy = 0 

-b/2 -a/2 y 


/ b/2 / a/2 [ N( B ) - J(w, ) + L(w, )] e dxdy = 0 

-b/2 -a/2 y 

From integrating Eqs. (71) and (72), one can find that the constants 
and B 2 of the slope functions are found to be 



-11* b 10 (i) 2 . (b 5 .b 12 ) ( I) 2 ] 
a a b 

„ 4 _ 2 -2 -4 -2 ,-,2 

i - k, (i) -k 3 fi) (i) - k 5 fi) + k 6 (*) + k e a 

a abb a b 


b 2 = 


— [1 + b 3 fl) . (b, . bj) (-) j 
b b a 


ki (-)*- k 3 ft* ft' - k 5 ft 4 ♦ k 6 (-/ * k 8 (l) 2 
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Note that the constants Bi and B 2 depend on plate dimensions and 
coefficients b. and k. . With no transverse shear effects, T c = ( 


coefficients b. 




and k^ vanish and the slope functions reduce to 
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It is obvious that the expressions for B 2 and B 2 given in Eqs. (73) and 
(74) are consistent with no transverse shear condition given in Eqs. (75) 
and (76) . 

Stress and Strain Response 

The strains at any point of the laminate are given by equation (2) and 
by using the constitutive relationships given in Eq. (10). One can write 
Eq. (2) as 
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(77) 


where z is thickness coordinate. Once the strains are known the stresses 
in kth layer can be determined from 
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where Q.^ are the transformed reduced stiffnesses for kth layer. 
Suostituting for stress function F as given in Eq. (44) and slope 
functions a and 6 as given in Eqs. (42) and (43), a general expression 
for strains at any point in the structure is given as 
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C x = An*n x + Ai 2 *n y + Ai6* n X y 

C 2 y = Al2*n x + A22*riy + A26* n X y (82) 

C 2 xy = A 16 *n x + A26*riy + A66*^ X y 

with 
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The mean-square strain is then related to mean-square modal amplitude in a 
general expression as 


E [ e 2 ] = C? E[ q 2 ] + 3 & (E[ q 2 ]) 2 (84) 


where E[ q 2 ] is mean-square displacement of Eq. (70). By making use of Eq. 
(78), an expression similar to Eq . (79) can be obtained for stress at any 
point in the structure, as given by 
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and 



It is obvious that a similar expression to Eq. (84) can be employed for 
computing mean-square stresses. 

Transverse Shear Stresses 

The transverse shear stresses can be obtained either by using the con- 
stitutive equations f (Eqs. (2) and (5)) or by integrating equil ibri un equa- 
tions (of three-dimensional elasticity in the absence of body forces) with 
respect to the thickness coordinate. Reddy (Ref. 33) felt that the second 
approach not only gives single-valued shear stresses at the interfaces but 
yields excellent results for all theories in comparison with the three- 
dimensional solutions. In view of its accuracy, inspite of the fact that 
the use of stress equiliDrium equations in the analysis of laminated plates 
is quite cumbersome, the integration of equilibria equations is used for 
the determination of shear stresses t X 2 and Xy Z< Thus, integrating the 
equilibrium equations with respect to the thickness coordinate z yields 

T xz = -\ n fo x,x + T xy,P dz + < 88) 
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Using Eq. (77) and the inplane stress-strain relationships given in Eq. 
(78), Eqs. (88) and (89) can be written as 
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(k) 

where Q. . = Q.. ' but for brevity the superscript k is omitted in Eqs. 

J 1 J 

(90) and (91) and hereafter. Expressions for n y , n and n yi/ are given 

x y xy 

in Eq. (83) and are defined in Eq. (81). Equations (90) and (91) are 
integrated with respect to the thickness coordinates z and rearranged to 
give 


T xz = D lX z 9 + 0 2XZ q 2 + f(x,y) 


(92) 


T yz = Dl yz q + ° 2 yz q2 + g ( x »y) 


(93) 


31 


where the functions f(x,y) and g(x,y) have to be determined from 
continuity considerations and 
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Equations (92) and (93) provide general expressions for the transverse shear 
stresses at any location (x,y,z). Since the excitation is random and the 
interest is in determining the mean-square x xz and x . Before the mean- 
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square stresses are calculated, the functions f(x,y) and g(x,y) are to 
be evaluated. 

The following scheme explains the procedure for evaluating the func- 
tions f(x,y) and g(x,y) for a four layered symmetrical laminate with the 
stacking sequence of (0/90/90/0), see Fig. 2. Because of the symmmetry 
only two layers are considered. Numbers 1, 2, 3 refer to the thickness 
coordinate z. Superscripts I, II refer to layer numbers as well as lamina 
material properties. Considering layer I and the point z (1), one can get 
from Eq. (92) 


xz( I) = D lxz(l) q + D 2xz( 1) q2 + f ( x ’ y) 


( 101 ) 


Now imposing the condition that T xz = 0 on top and bottom of the laminate 
results in 


fI < x ,y)= -D lx ‘ (1) q - l) 2x ‘ (1) q* 


( 102 ) 


By considering the point 2 of layer I (which is an interface point for laye 

I and layer II), from Eq . (92) one gets 


xz(2) ‘ U lxz(2) q U 2xz( 2) 
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The substitution of Eq. (102) in Eq. (103) results in 
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The condition that x xz = 0 at point 1 can be expressed as 
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1 r 

T xz(l) " b lxz(l) 


(105) 


* + G 2XZ(1) 


q2 


with G ixz(l) = G 2 xz(2) = 0. 

Now Eq. (104) can be written as 

t xz(2) = G 1xz( 2) q + G 2 xz{ 2) q2 

where 

G 1xz(2) = D 1xz(2) ' °lx I z(l) 

G 2xz(2) = °2xz(2) ' °2xz(l) 

Now considering the same point 2, but in layer II one can write 

t xz(2) = °lxz(2) q + D 2xz( 2) q2 + 

since continuity consideration requires this 

II I 

t xz(2) " t xz(2) 

using Eqs. (106), (109) and (110), one can show that 
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Next considering point 3 in layer II, one gets 


11 n II ^ n 11 2 JI V 

xz(3) = D lxz(3) « * D 2xz( 3) ’ + f ' x ^> 


( 112 ) 


substituting for f^fx.y), Eq. (112) becomes 
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t xz(3) 


^°1xz(3) ' D 1xz(2) + G 1xz( 2)^ q 
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Equation (113) can be rewritten as 


II 


T xz(3) = G 1 xz( 3) q + G 2 xz( 3) q 


(114) 


Once the constants Gixz, G 2 xz are evaluated at points 1, 2, 3 etc., along 
thickness, the mean-square t xz can be evaluated using an expression 
similar to Eq . (84). In similar manner the functions g(x,y) is evaluated 
and the constants Giy Z » G 2 y Z are determined for calculating mean-square 
T yz- 

NUMERICAL RESULTS AND DISCUSSION 

Numerical results for a symmetric cross-ply plate are presented first. 
A four-ply square laminate (12x12 in.) with layers of equal thickness and 
subjected to a uniform random pressure is considered. The plate is simply 
supported on all four edges. For the examples presented, a representative 
high-modulus graphite-epoxy with the following material properties is used. 


Ei/E 2 = 40 
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G23/E2 = 0.5 


G12/E2 = Gi 3 /E 2 = 0.6 

v 12 = 0.25 

E 2 = 0.75xl0 6 psi 

p = 2.4xl0- 4 lb-sec 2 /in . 4 


Table 1 shows the nondimensional fundamental frequency u> o = 

(u 0 a 2 /h) yjp/E 2 as a function of plate length to thickness ratio (a/h). 

The linear small deflection plate theory with and without transverse shear 

deformation is used. The Navier series solutions obtained by Reddy and Phan 

(Ref. 41) are also given in Table 1 for comparison. It clearly indicates 

that the present method gives good frequency predictions. 

Table 2 shows the RMS nondimensional maximum deflection W , and RMS 

max 

nondimensional maximum stress in the major material direction "aj , versus 
plate length to thickness ratio for the same cross-ply laminate at 130 dB 
(Ref. 2xl0- 5 N/m 2 ) sound spectrum level (SSL). Stress cFj is measured at 
(0,0, h/2). Table 3 shows the nondimensional equivalent linear frequency 
IT = (na 2 /h) / p / E 2 versus plate length to thickness ratio. Examination of 


Tables 2 and 3 reveals that for moderately thick plates (a/h<20) the small 
deflection theory with shear deformation (x = 0 and T $ = 1) and for thin 
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plates (a/h>50) the large deflection theory without shear deformation (x * 0 
and T $ = 0) would give accurate predictions as indicated by agreement with 
the theory including both nonlinear and shear effects (X * 0 and T $ = 1) on 
maximum deflections, inplane stresses and equivalent linear or nonlinear 
frequencies. This is clearly evident from Figs. 3 to 5, where the values 
given in Tables 2 and 3 for RMS "w , RMS stresses and frequency TT are 
shown plotted against a/h. 

The RMS nondimensional transverse shear stresses ai* and F 5 are de- 
termined by integrating three dimensional equilibrium equations as explained 
earlier, a" 4 is measured at (0,b/2,0) and F 5 is measured at (a/2,0,0). 
Table 4 shows ~a h and cT 5 values calculated for different a/h ratios at 
130 dB SSL. Figure 6 shows the transverse shear stresses 7j 4 and "ct 5 across 
the plate thickness for the laminate under consideration with a/h = 10 and 
at 130 dB sound spectrum level. It is evident from Fig. 6, that for thick 
laminates, the plate theory without transverse shear over predicts the RMS 
transverse shear stresses compared to plate with shear deformation. 

Figures 7, 8 and 9 show the RMS (maximum deflection/h) , RMS maximum 


m a v-r-n a I 


.tress cfj and equivalent linear frequency n, respectively, for 


the simply supported four-layer cross-ply square plate with a/h = 200 at 
sound spectrum level varying from 90 to 130 dB (ref. 2xl0- 5 N/m 2 ). Results 
shown are using the four formulations discussed earlier and there is no 
appreciable difference between the results using theories with and without 
transverse shear. The linear and nonlinear solutions agree at low values of 
SSL, but disagree at high values. For acoustic excitations of sound 
spectrum levels less than 90 dB, the small deflection assumption will give 


good predictions for the composite panels studied as the linear and 
nonlinear solutions coalesce. At 130 dB SSL, however, the small deflection 
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theory predicts that the plate would deflect to a value of 6.6 times of 
plate thickness; whereas, the large deflection theory (with or without 
shear) gives the much smaller value of 1.6h. This smaller value seems more 
reasonable, intuitively. Similarly, for the RMS stresses and equivalent 
linear frequency, at high SSL values (for thin laminates) the small 
deflection theory over predicts the RMS deflection and stresses, and under 
predicts the frequency as compared with large deflection plate theory. 

Finally, the nondimensional ized fundamental frequencies ai 0 for a 
symmetrical angle-ply (6=±45) square plate (12xl2in.) for different a/h 
ratios and number of layers are shown in Table 5. The material is graphite- 
epoxy and the material properties are those that were given in Ref. 42. The 
nondimensional ized frequencies compare very well with the values given in 
Ref. 42. 


SUMMARY AND CONCLUSIONS 

The main objective of this study is to predict mean-square inplane 
stresses and transverse shear stresses that develop in symmetrical composite 
laminates when they are subjected to acoustic excitation. For moderately 
thick laminates, the transverse shear effects are considerable. In this 
report, equations of motion are developed which include geometric large 
amplitude nonlinear effects (von Karman theory). The transverse shear de- 
formation effects are included. By various operations, the slope functions 
are eliminated from the equations of motion, which are expressed in terms of 
stress function F and displacement function w. These equations can be 
considered as an extension of von Karman' s nonlinear equations of plates. 
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In view of the complexity of the equations, a single-mode Galerkin 
procedure is employed to obtain a nonlinear modal amplitude equation for the 
forced vibration of plate. The excitation is assunmed to be stationary, 
ergodic and Gaussian with zero-mean. The equivalent linearization method is 
employed. The fundamental frequencies that are obtained for both cross-ply 
and angle-ply laminates found to be in good agreement with those that are 
available in the literature. Root-mean-square deflections and RMS inplane 
stresses are calculated. Using the three-dimensional equilibrium equations 
and continuity considerations, the RMS transverse shearing stresses in the 
laminate are determined. 

The effects of transverse shear are considerable if the plate lengths 
are less than 20 times the thickness (a/h<20). These effects should not be 
ignored in moderately thick and thick plates and small-deflection theory 
with shear deformation would give accurate predictions of maximum 
deflection, frequency and inplane stresses. For thin plates (a/h>50), the 
linear and nonlinear solutions agree at low values of SSL, but disagree at 
high values. The small deflection theory over predicts the RMS deflection 
and stresses, and under predicts the frequency as compared with large 
deflection plate theory. The large deflection theory with transverse shear 
effects neglected would give accurate predictions for thin panels at high 
SSL values. For a particular value of a/h, therefore, one of the three 
simpler theories can be chosen that provides accuracy equal to the more 
cumbersome theory that includes Doth shear and large deflection effects. 

The prediction of transverse shearing stresses is required for 
understanding sonic fatigues of composite laminates, especially the unique 
inter-laminar failures. 
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Table 1. 

Nondimensional 
Supported Four 

Fundamental Frequency (&» 0 a 2 /h) 
-Layer Cross-Ply Square Plate. 

Jp/E 2 of a Simply 

a 

No Shear 

Shear 

h 

Navier 41 

Present 

Navier 4 1 

Present 


Solution 

Result 

Solution 

Result 

5 

18.215 

18.891 

10.820 

11.554 

10 

18.652 

18.891 

15.083 

15.662 

20 

18.767 

18.891 

17.583 

17.872 

50 

18.799 

18.891 

18.590 

18.715 

100 

18.804 

18.891 

18.751 

18.846 

200 

- 

18.891 

- 

18.880 

400 

- 

18.891 

- 

18.888 

1000 

- 

18.891 

- 

18.890 
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Table 2. 


Nondimensional RMS Maximun Deflection 10( RMS w)E 2 h 3 / (a 4 \J(PSD)f n ) 

and RMS Maximum Stress (RMS 01 ) h 2 /(10a 2 V ( PSD) f«) of a Simply 
Supported Four-Layer Cross-Ply Square Plate at 130 dB Sound 
Spectrum Level . 


a 

h 

Small Deflection 
No Shear Shear 

Large Deflection 
No Shear Shear 


RMS 

Maximun Deflection 


5 

0.4026 

0.2471 

0.4026 

0.2471 

10 

0.4026 

0.3478 

0.4026 

0.3478 

20 

0.4026 

0.3873 

0.4026 

0.3873 

50 

0.4026 

0.4001 

0.3993 

0.3968 

100 

0.4026 

0.4019 

0.2810 

0.2806 

200 

0.4026 

0.4024 

0.0973 

0.0972 

400 

0.4026 

0.4025 

0.0293 

0.0293 

1000 

0.4026 

0.4026 

0.0059 

0.0059 


RMS 

. Maximun Stress 



5 

0.8009 

0.1521 

0.8009 

0.1521 

10 

0.8009 

0.4467 

0.8009 

0.4467 

20 

0.8009 

0.6780 

0.8008 

0.6780 

50 

0.8009 

0.7789 

0.7951 

0.7732 

100 

0.8009 

0.7953 

0.5923 

0.5885 

200 

0.8009 

0.7995 

0.3284 

0.3280 

400 

0.8009 

0.8005 

0.2788 

0.2788 

1000 

0.8009 

0.8008 

0.2739 

0.2739 
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Table 3. Nondimensional Equivalent Linear Frequency (fi a 2 /h)Vp/Ej of a 

Simply Supported Four-Layer Cross-Ply Square Plate at 130 dB Sound 
Spectrin Level. 


a 

h 

Small 
No Shear 

Deflection 

Shear 

Large 
No Shear 

Deflection 

Shear 

5 

18.891 

11.554 

18.891 

11.554 

10 

18.891 

15.662 

18.891 

15.662 

20 

18.891 

17.872 

18.900 

17.872 

50 

18.891 

18.715 

19.047 

18.870 

100 

18.891 

18.846 

27.063 

26.995 

200 

18.891 

18.880 

78.196 

78.146 

400 

18.891 

18.888 

259.463 

259.396 

1000 

18.891 

18.890 

1287.297 

1287.264 
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Table 4. Nondimensional RMS Maximun Transverse Shear Stress 

(RMS t )h/10a V7PSD] f ) and (RMS t )h/10a V (PSD) f n ) 
of a Simply Supported Four-Layer Cross-Ply Square Plate°at 
Sound Spectrun Level . 


130 dB 


a/h 

Small 
No Shear 

Deflection 

Shear 


Large Deflection 
No Shear Shear 

5 

0.1906 

RMS Maximum 
0.0797 

(<*<♦) 

0.1906 

0.0797 

10 

0.1906 

0.1459 


0.1906 

0.1459 

20 

0.1906 

0.1774 


0.1906 

0.1774 

50 

0.1906 

0.1883 


0.1890 

0.1868 

100 

0.1906 

0.1900 


0.1330 

0.1327 

200 

0.1906 

0.1905 


0.0461 

0.0460 

400 

0.1906 

0.1906 


0.0139 

0.0139 

1000 

0.1906 

0.1906 


0.0028 

0.0028 

5 

0.4954 

RMS Maximum t xz 
0.0966 

0*5 ) 

0.4954 

0.0966 

10 

0.4954 

0.2786 


0.4954 

0.2786 

20 

0.4954 

0.4203 


0.4953 

0.4203 

50 

0.4954 

0.4819 


0.4913 

0.4779 

100 

0.4954 

0.4919 


0.3458 

0.3434 

200 

0.4954 

0.4945 


0.1197 

0.1195 

400 

0.4954 

0.4951 


0.0361 

0.0361 

1000 

0.4954 

0.4953 


0.0073 

0.0073 
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Table 5. Nondimensional Fundamental Frequency (co 0 a 2 /h)Vp/E 2 of a Simply 
Supported n-Layer angle-ply (6 =±45) square plate. 


a/h No. of Layers 

No Shear 

Ref. 42 Present Result 

Ref. 42 

Shear 

Present Result 

5 

3 

25.82 

25.82 

12.78 

12.21 

10 

3 

25.82 

25.82 

19.38 

18.89 

10 

5 

25.82 

25.82 

19.23 

19.06 

10 

7 

25.82 

25.82 

19.19 

19.12 

20 

3 

25.82 

25.82 

23.62 

23.40 
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Fig. 3. Nondimensional RMS maximum deflection of a simply supported 
four-layer cross-ply square plate at 130 dB sound spectrum 
level. 
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dB sound soectrum level. 



Linear, no shear 
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Frequency in Cycles per Second, a/h 

. Nondimensional equivalent linear frequency of 
simply supported four— layer cross — ply square 
plate at 130 dB sound spectrum level. 



z/h z/h 


0.5 


0.3 

0.1 

- 0.1 

-0.3 

-0.5 

0.5 

0.3 

0.1 

- 0.1 

-0.3 

-0.5 



Fig. 6. Nondimensional RMS transverse shearing 

stresses of a simply supported four — layer 
cross-ply of a /h = 10 at 130 dB sound 
spectrum level. 
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Sound Spectrum Level, dB 
Re 2 x 10” 5 N/m 2 

Fig. 7. RMS maximum deflection versus sound spectrum 
level for a simply supported four — layer cross- 
ply square plate with a/h = 200. 



RMS Max. Stress 


Sound Spectrum Level, dB 
Re 2x10~ 5 N/m : 

Fig. 8. RMS maximum stress versus sound spectrum 
level for a simply supported foui — layer 
cross — ply square plate with a/h = 200. 





Sound Spectrum Level, dB 
Re 2x10“ 5 N/m 2 

Fig. 9. Equivalent linear frequency versus sound 

spectrum level for a simply supported four- 
layer cross— ply square plate with a/h = 200. 
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appendix 


The boundary conditions for clamped plate are 
x = ± a/2 : w = a = g = 0 


(Al) 


y = ± b/2 : w = a = 8 = 0 

The deflection function which satisfy the boundary conditions on all 
edges is asstmed as 

q(t)h , 2 ttx Wi , 2 iiyx 

w = ( 1 + cos J ( 1 + cos i) 

4 a b 

The stress function is of the following form 

F = F + F 
c p 

in which the particular solution is 

q 2 k2 r 2 

F = - — — [C 10 cos X + C 01 cos Y 

p 32 

+ Cn cos X cos Y + C 20 cos 2 X 

+ C 02 cos 2 Y + C 21 cos 2 X cos Y 

+ C j 2 cos X cos 2 Y + S u sin X sin Y 


four 


(A2) 


(A3) 
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(A4) 


+ S 21 sin 2 X sin Y 
+ Si 2 sin X sin 2 y] 


where 


X = 


2*x 

a 



a/b 


(A5) 


and the constants C.. and S. . can be expressed in terms of the inverted 

' J * J 

extensional laminate stiffness and length-to-width ratio, r, of the panel 
as 


ClO - 1/ A*22 
Coi = l/(r 4 A*n) 

Cn = 2Gs/(G5 - Gg) 
C 2 o = 1/(16 A* 22 ) 

C 02 '= 1/(16 A*i 1 r 4 ) 

, 2 2 ^ 

Ci 2 = 63/(63 " G4) 

2 2 

C21 - Gi/(Gi - G2) 
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2 2 

Si i = - 2Gs/C Gs - G6 ) 

2 2 

Si 2 = - Gi^/C G 3 " G4 ) 

2 2 

S 21 = - G 2 /( Gi - G 2 ) 

Gi =16 A *22 + 4(2A*i2 + A*66 )r 2 + A*ii r 4 

G 2 = 16 A *26 r + 4 A*i 6 r 3 

G 3 = A *22 + 4 (2A*i2 + A*66 ) r 2 + 16 A*n r 4 

G 4 = 4 A *26 r + 16 A*i 6 r 3 

G 5 = A *22 + ( 2 A*i 2 + A* 66 ) r 2 + A*n r 4 

G6 = 2A*26 r + 2A*i6 r 3 (A6) 

The complementary solution is of the sane form as given by Eq. (49) with 

# = 3q 2 h 2 tt 2 / A 22 * - Ai 2 * \ 

32(A 11 *A 22 *.Ai2* 2 ) \ a 2 b 2 / 

(A7) 
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3q 2 h 2 ir 5 


An* . A 12 *^ 


y 32( An* A 22*" A 12* 2 ^ \ b2 a2 


N xy = ° 


The modal amplitude equation is of the same form as that of Eq. (52) with 


-5 - - {3D 11 (i) , » + 2(D 12 + 2D 66 )(!)2 (1)2 + 3D 22 (1) 4 
m a abb 


- 4 [3Pi(l)6 + P 3 (!L)4 (1)2 + r 5 (!)2 (1)4 + 3P 7 C— ) 6 ] } (A8) 

a a b a b b 


and the modal mass m is given by 


m = 5£ii {1 - i + i k, (i) 2 (l) 2 * 4k s (i)* - k 6 (I)4- k s (l)2]) ( A9) 

16 3a3ab b a b 


The nonlinearity coefficient is of the same form as given in Eq. (54) with 


X D * ^ f C l° f 1 - 8k l r1 ) 4 - 8k 3 f-) 2 (^ 2 - 16k 5 (— ) 4 + 2k 6 (— ) 2 * 4k e (— ) 2 ] 
p 16mb 4 a a b b a b 


2 . 2 


+ C 0 ifl-16ki(l) - 8k 3 (-) (*-) - 8k 5 (— ) + 4k 6 (l) + 2k 8 (L) ] 
a a b b a b 


2 „ 2 


t Ci ,f l-8k 1 f— ) - 8k 5 (— ) fl) + 2k 6 (— ) + 2k 8 (— ) ] 
a a b a o 
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+(C 20 aC 02 ) [1-16 k^-)* - 16k, (1)2(1)* . 16k 5 (*)- a 4k 6 (")2 a 4kg(— ) 2 ] 

a a b a a b 


tIc 21 [l-16k,(L)4 - 8k, (—)* (— )* - 8k s (l)k + 4k 6 (I)2 + 2k 8 (i)2)] 
2 a a b b a b 


* i Cj 2 [l-Sk 1 (— )■*- 8k 3 (l)2(!L)2 - 16k 5 (— ) 4 a 2k 6 (£)2 a 4k e (I)2)] 
2 a a b b a b 


* <S 21 aS 12 )[4k 2 (L)3(l) a 4k„(’)(i)2 - k y (_)(_)]} 
a b a b a b 


(A10) 


A = f (A 22 * - A 12 *r2) (l-4f4k jfi) 1 * 

c 128 m (A n *A 22 *-A 12 *2)b4r4 a 


+ i k 3 (!)2(L)2 + 1 k 5 (l)- . k 6 fl)2 . 1 k 8 fI)2])+(A n *r4-A 12 *r2) 
3 a b 3b a 3b 


(l-4[i k,(!) 4 a i k, (— )2 (— )2 a 4k, (1)* - 1 k 6 (!) 2 - k,(l)2])J 
3a3ab b3a b 

with transverse shear effects neglected, ie T s = 0. The coefficients 
P. are zero, Eqs. (A8)-(A11) reduce to 

“o = t 3 + 2(D i2 + 2D 66 ) r2 + 3 D 22 r4] 

9pha H 

9ph 2 
m = 

16 

U I 

^ - — fCjo + Cqi + Cji + C 02 + C 20 + — (C 2i + Cj 2 )i 

K 9p b H 2 


(All) 



(A12) 


(A13) 


(A14) 
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X = lUil ( A 22 * ~ 2 A 12 * r 2 + An* r** 

C 8 p a 4 \ A 1X * A 22 * - A 12 * r 2 

Equations (A12) - (A15) are the same as in Ref. 3. The slope functions a 
and $ for the clamped case are assumed to be 


« (x,y,t) 

= Bj q(t) sin — - (1 + cos ) 

a b 

(A16) 

3 (x,y,t) 

= B 2 q(t) (1 + cos 2Tr ^) sin 2ir,y - 

a b 

(A17) 


with 


ill [1 t 4bi„ (1)2- i (bs- b 12 ) (1)2] 

2 a a 3 b 

Bj = (A18) 

U - 4 f 4k, (1)4 ♦ i k,(!)2(l)2 * i kg (1)4 - k 6 (lf - 1 k,(l)2]] 
a3ab3b a3b 

— Tl + 4 b 3 (I) 2 - i (bs-bj) (1)2] 

2b V 3 a 

B 2 = (A19) 

fl - 4 ri kjfLf + 1 k 3 (— ) 2 — ) 2 + 4k 5 (L)-- 1 k 6 (2L)2 - k 8 (!)2]} 

3a 3ab b3a b' 

The S i of Eq. (80) are defined as 

s, * — - B, 
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( A20) 


2tt 

s 2 = — B 2 

b 


S 3 = -2(1 Bj + 1 B 2 ) 
b a 


The of Eq. (82) are defined as 


= i (!HL) 2 (4S 12 sin X sin 2 V + S 2 i sin 2X sin Y + S n sin X sin Y 


8 b 


+ 4Ci 2 cos X cos 2 Y + C 2 i cos 2 X cos Y + 4Co2 cos 2 Y 


+ Ci i cos X cos Y + C 0 1 cos Y) + 


3h 2 n 2 


A22* Al 2* ' 


32(An*A 22 * - Ai 2 * 2 ) \a 2 b 2 


n =1 (JZlZL) 2 (S 12 sin X sin 2 Y + 4S 21 sin 2 X sin Y + S n sin X sin 
y 8 a 


+ Ci 2 cosX cos2Y + 4 C 2 i cos2X cos Y + 4 C 2 o cos2X + Ci i cosX cosY 


+ C 10 cosX) + i-1- 2 


An* Ai 2 * 


32(Ai i*A 22 *-Ai 2 * 2 ) \ b 2 a 2 


Sab 


L_ (hnr) 2 (2S 12 cos X cos 2 Y + 2S 21 sin 2X cos Y 


(A21) 

Y 


(A22) 
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+ Sji cos X cos Y + 2Ci2 sin X sin 2Y + Z 2l sin 2X sin Y 


+ Ci i sin X sin y) 

The analagons equations for Eqs. (94) - (97) are 

° lX z = °- 5z2 [( f i" f 2 ) sin X cos Y + (F 3 .F 4 ) cos X sin Y + E 10 sin X 
+ E 30 sin y] 

D 2xz = - z[Ri Qn + «2 Q12 + ( R 3 + M Qie + R 5 Q 26 + R 6 Qee] 

D 1 yz = °* 5Z2 [ ( F s ) sin X cos Y + (F 7 -F 8 ) cos X sin Y + E 50 sin X 
+ E70 sin y] 

y Z = - Z f R l Ql6 + (R 2 + R fi ) Q 26 + R 3 Q66 + R 4 Qi 2 + R 5 Q 22 ] 
where 

F i = — HTu Sj + I 7 12 S 2 ) 

d 

F 2 = ~ (^66 S 3 ) 

D 

F 3 = — (^16 $1 + 7 2 6 S 2 ) 


(A23) 

( A24) 
(A25) 

(A26) 

(A27) 
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(A28) 


F 4 = fQi 6 S 3) 

a 

F 5 = — f$16 S 1 + ^26 S 2 ) 
3 

^6 = — (^26 S 3 ) 
b 

f 7 = ?1 (TJ 12 Sj + U 22 s 2) 

b 

F 8 - i- (^66 S 3 ) 

a ' 

E 1 0 = f^ii s i) 

3 

E 30 = 0?26 S 2 ) 

b 

E50 = — (Qie Si) 

a 

E 70 = (If 22 S 2 ) 

b 
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